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JrJ ■ Abstract 

(—1 

■ ■ We consider the interference of two overlapping ideal Bose-Einstein condensates. The usual 

description of this phenomenon involves the introduction of a so-called condensate wave functions 

having a definite phase. We investigate the origin of this phase and the theoretical basis of treating 

^5 [ interference. It is possible to construct a phase state, for which the particle number is uncertain, but 

Q ' phase is known. However, how one would prepare such a state before an experiment is not obvious. 

We show that a phase can also arise from experiments using condensates in Fock states, that is, 

having known particle numbers. Analysis of measurements in such states also gives us a prescription 

^ ■ for preparing phase states. The connection of this procedure to questions of "spontaneously broken 

f"^ . gauge symmetry" and to "hidden variables" is mentioned. 

o 
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I. INTRODUCTION 

One of the most impressive experiments using trapped Bose gases is the interference 
experiment of Ketterle and co-workers.— Two condensates are are separately prepared and 
allowed to overlap. An interference pattern then arises showing the remarkable quantum 
coherence of the condensates. There have been other interesting condensate interference 
experiments as well.—"— If one assumes that the separated clouds initially have a definite 
phase relation, then the experiments are well-described by straight-forward theory.^ However, 
questions immediately arise. Do separately prepared condensates have a phase relation?— 
The preparation of the sample certainly did not involve the establishment of a state with 
known phase. More likely the particle number in each cloud would have, or could have, 
been initially known ahead of time. Nevertheless an interference pattern emerges with some 
well-established phase. So how is this possible? The question was answered in several 
theoretical papers showing how a phase appears even when the two clouds are prepared in 
Fock states, that is, states with sharply known particle numbers.—"— In this paper we revisit 
the question and present a derivation of this result. This result is quite satisfying, since it 
in a sense justifies the usual simple assumption of the interf erring coherent systems having 
a well-defined, but completely unknown, initial phase relationship. 

In order to discuss the properties of condensates and superfluids one usually introduces 
theoretically a so-called "order parameter" or "condensate wave function" defined as ( V'(r) ) 
where '?/'(r) is an operator destroying a particle at position r. The resulting quantity has a 
magnitude equal to the square root of the condensate density and also a phase. Obviously 
the state in which (ip{r)) is non-zero cannot have a fixed particle number. With one 
of these wave functions for each condensate, it is straightforward to discuss interference 
of the two, since each one has its own phase, and an interference pattern arises with a 
relative phase equal to the difference between the individual phases. In essence one has 
described each condensate by a single-particle wave function so that interference is no more 
than the overlap and interference of two classical waves. However, one can ask how this 
single-particle wave function might have arisen. Indeed its existence involves the question 
of "spontaneously broken gauge symmetry"—"—, the necessity of which has been brought 
into question in recent years.—"— 

Suppose we consider a condensate described by a wave function e'^^^'^'^'^\ One might de- 



scribe the direction specified by the angle by a "spin" in a two-dimensional plane. How 
do we prepare such a state? What is it that picks out the direction of this pseudo-spin from 
the all the degenerate possible directions? There is an analogy with ferromagnetism, where 
there is symmetry in the possible degenerate three-dimensional directions of the magneti- 
zation. Somehow, say, because of a small external field, a particuliar direction is space is 
picked out for the magnetization. In a similar way the phase angle is picked out. One can 
suppose in the ferromagnetism case that, in actual practice, there is always some small field 
to pick out a preferred direction, so that the symmetry is broken. However, the local field 
that is used theoretically to pick a phase direction does not exist in nature. Indeed one 
would have to prepare a state that does not have a fixed particle number. The treatment 
of phase (actually relative phase) in Sec. IV, and the spontaneous appearance of a relative 
phase under the effect of measurement of particle position in Fock states, avoids violating 
particle conservation and does not require use of any symmetry-breaking field and so helps 
in this regard. 

A closely related feature of the discussion is the idea that the phase emerging from 
successive measurements of particle position starting with a Fock state is somewhat like 
the emergence of a "hidden" or "additional" variable in quantum mechanics.— »^>^ Was 
the phase there before the experiment started, or did the experiment itself cause it to take 
on its final value? Hidden variables can be invoked to specify non-commuting variables. 
In standard quantum mechanics, particle number and relative phase can be considered 
conjugate variables; the knowledge of one excludes that of the other. As we measure particle 
position we will see that the knowledge of the particle number becomes less certain while 
the uncertainty of relative phase decreases. 

In the next sections we first discuss the kind of state that has a known phase. Obviously 
with this state the interference pattern emerges with just the prepared phase. Among 
these states are are the coherent states of Glauber.i^ These can either have particle number 
completely unknown or have the total number of particles in two clouds known (in which 
case they are called "phase states"), although the number in each cloud is still unknown. 
We then derive the interference pattern starting with Fock states and see the emergence of 
a relative phase even though no phase was present at the beginning of the experiment (or 
was at least "hidden"). We even find a way to prepare a state that has a known relative 
phase. The controversal theoretical constructs are seen to be unnecessary. 



II. SIMPLE VIEW OF AN INTERFERENCE PATTERN 

A gas (or liquid) undergoing Bose-Einstein condensation (BEC) is often described by 
a classical field known as an "order parameter" or "condensate wave function." Such a 
quantity can arise in several ways. Suppose that ip{r) represents a second-quantized operator 
destroying a boson at position r. Then the one-particle density matrix is defined as p(r, r') = 
ijj (r)^(r) ) . Penrose and Onsager^^ showed that a criterion for a Bose condensate or "off- 



diagonal long-range order" is that the density matrix have the form 

p(r,r') = 'r(rMr') + /(r,r') (1) 

where /(r, r') vanishes when |r — r'| ^ oo . The function 4^{r) is the condensate wave 
function. One often assumes the system in a state such that the destruction operator itself 
has a non-zero average: 

(^(r))=^(r) = yV^e^^W (2) 

where no is the condensate density and its phase. Such a state is said to have "spon- 
taneously broken gauge symmetry" because a particular phase (out of many possible de- 
generate phase states) has been chosen.—"— The Gross-Pitaevskii equation is a non-linear 
Schrodinger equation for ip{r), which has been used extensively, with remarkable success, to 
describe interacting trapped BEC gases within mean-field approximation.— 

In describing the interference pattern of the experiment of Ref. [l| one has to consider 
the overlap of Bose clouds released from harmonic oscillator traps.— This leads to some 
interesting features such as fringes whose separation changes with time. In our analysis here 
we will consider only plane waves and ignore any time evolution. Thus suppose we have an 
order parameter that involves two condensate clouds, having condensate densities n^ and rib 
in momentum states k^ and k^. This dual order parameter has the form 

^(r) = ^e'^-'e''^'^ + y/^te'^'^'' e'^" . (3) 

The density of the combined system is then 

n(r) = I y/^e'^'^^'e'^'^ + ^n^e^'^^'e^'^' f = n{l + x cos(k ■ r + 0)) (4) 

where n = n^ -|- n^, k = k^ — kf,, = 0^ — 0^, and x = 2-^fn^jiiln. We have an interference 
pattern with relative phase k ■ r-|-0. The phase shift is measureable although the individual 
phases 0„ and 0^ are not. 



This analysis is simple, but it requires the preparation of the system in a state with 
known individual phases. How do we do that? What is the nature of such a state? Clearly 
the expectation value of Eq. Q cannot be in a state of definite particle number or the 
expectation value would vanish. Next we investigate this question more deeply. 

III. PHASE STATES 

As noted by Johnston^i in this journal, the coherent states introduced by Glaubei^ for 
photons are appropriate for superfluids. These well-know states have been reviewed in this 
journal on occasion, ^^ and appear in texta^S as well. They are also called "classical states" 
and are the minimum uncertainty states of the harmonic oscillator.— Here we do not use 
them in full generality, but rather use a subset of them known as "phase states." (Those 
wishing to see the full treatment of coherent states in a treatment of a condensate wave 
function are referred to Appendix A.) Phase states describe two condensates (in states k^ 
and kfe) with variable particle numbers, Na, Nt (both macroscopic), but fixed total number 
N = Na + Nf,. No other momentum states are occupied. If particle creation operators a^ 
and 6^ (obeying Bose commutation relations) for the two states act on the vacuum to put 
particles into these two states, then we define the (properly normalized) state as 
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where we define the quantities «» as complex and separate them into magnitudes 7^ and 
phases 0j according to the notation 

«i=7^e''^^ (6) 



Also g = \/^\ + 7^- We can easily compute the average number of particles Na in this state. 
Use the first form of Eq. (jSJ to give 



a|a,a,;iV) = ^y \ ^' ^^ ^, a^-a^-^" v^ |iV, - 1, iV - iV,) 
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where we have taken A^^ = A^^ — 1. Thus 

Na = {aaab;N\a^a\aaab;N) =^1^^ — 2 (8) 

Ta + lb 

Similarly we get Nb = llN/i^l + 7^) so that 7^ = \ai\ = \/Wi and 51 = A'^ + A'f, = A^- 

The fact that A^ is known in our phase state does not affect the results for interference 
patterns, which depend just on relative phase. Indeed such states have been used on many 
occasions to discuss the interference of two condensates.— »^»^>^ Our state can be used to 
discuss the how relative phase can be conjugate to particle number. Write it in a form that 
makes the phases explicit: 

V y {Na+N,=N) ViVa.iVfe. 

Now change the phases to relative phase = 0^ — 0^ and total phase $ = (1/2) (0^ + 0„) 
and then take the derivative with respect to : 

^ V ^ {Na+Ni=N) yi^a-^^b- 

= [a^a-b^b] \aaab; N) (10) 

The phase derivative operator gives the same result as the number difference operator so 
that and (A^ — Nb) are conjugate variables.—!^ Note that the total phase appears only as 
an external factor exp [iA^$] and so has no physical significance. Thus the individual phases 
have no physical significance; only the relative phase is a meaningful quantity. 

In order to emphasize this last point and put the phase state in a more compact form to 
treat interference we rename it and rewrite it as 

10, N) = -^^{a^ + je%Y |0) = -^= y W „ ,/ „ (76^^)^-" \n, N - n) (11) 
^Jg^N\ ^/g^ ~l\ ^K^ - '^V- I y J 



where now 7 = \/Na/Nb and (j) is the relative phase. Also now g = {1 + 7^)- We have simply 
dropped a meaningless factor of unit magnitude. 

Since we have just two occupied states, the terms in Eq. (jTrf) not referring to states ka 
and kb never contribute and we can more simply write 

V^(r) ^Cr=J^ (e^'^""a + e^'^'-'b) . (12) 

We can also make this more compact by writing 



yi (a + 6^*^-6) (13) 



with k = kfc — ka. The quantity extracted out e*'^"'" can again be dropped as a meaningless 
unit-magnitude factor. We want to consider how Cy. acts on the phase state. We have 

1 ^ / m 
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(■^ \n -l,N -n)+ e^'^'VA^ - n \n,N -n- l)\ (15) 



By changing variables in the first state to ra' = n — 1 we can put both terms in the same 
form so 

Cr|0,iV) = A(r,0)|0,iV-l) (16) 

where 



A(r,0) = J ^ (1 + Te^^e^*^-) (17) 

We can now easily compute the condensate wave function as 

(0, N\ ^(r) 10, N) = (0, N\ cr 10, N) = A(r,0) (18) 

If for the moment we put back the previously neglected leading phase factors we have 

(0, A^l V'(r) 10, N) = ^/W^e^^'^^'e''^'^ + Vfibe^^'^'e''^'^ (19) 

which has just the same form as Eq. Q, with average densities in place of precise densities 
and well-defined individual phases. However the individual phases are not measurable — only 
the relative phase is. Indeed the average density follows immediately as 

(0, A^l ^\r)V'(r) 10, N) = (0, N\ cjcr |0, A^) = |A(r, 0)|' = n (1 + a;cos(k ■ r + 0)) (20) 



just like Eq. (0]) where now n and x have obvious definitions in terms of averages. Thus a 
phase state provides a rigorous background for discussion of condensate wave functions and 
for the simphfied form of treating interference between the two condensates. How we might 
actually prepare one before the experiment is a separate difficult question, which we treat 
below. 

We will find it useful and necessary in the next section to consider more general cases 
in which we make measurements of many particle positions essentially simultaneously. This 
allows interference fringes to emerge where they would otherwise not occur. For our phase 
state such calculations are straightforward and add no additional information since the the 
multiparticle densities all factor in the phase states. For example, consider the expectation 
value of 

^ (r2)V'(r2)^ (ri)V^(ri) = cl^c^A^Cr^ ^ cl^cl^Cr^Cr,. (21) 

Of course, c^^ and cj don't commute, but in the approximation shown we are dropping a 
term of order A^ compared to one of order A^^. The last form is more convenient to use. By 
the eigenfunction behavior of the phase state we easily get 

{cP,N\i;\r2mvS\v^)i;iv^) |0,iV) ^ (0, iV| 4^4^Cr,Cr, |0,iV) = \A{r2)f |A(ri)|' 

2 

= Yin {1 + X cos{k ■ Ti + (j))) (22) 

This result generalizes to 

(<^, N\ ^\r„)^(r„) ■ ■ ■ ^Vi)^(ri) 10, N) ^ (0, N\ c^^ ■ • • 4^Cr„ ■ ■ ■ Cr, |0, N) 

m 

= JJn(l + xcos(k-r, + 0)) (23) 

Considering a state in the form Cr„ ■ ■ ■ c^ |\1/) is useful for interpreting an experiment. 
Here we can consider our experiment as detecting particle 1 and then 2 shortly thereafter, 
and so on. After m detections the wave function evolves to a state missing several particles. 
What is the nature of the state to which it has evolved? In the case of a phase state it is 
Cr„ ■ ■ ■ Cri 10, A^) ~ \'PiN — m) . However, it is more interesting to consider the case of |\1') a 
Fock state as we do in the next section. 



IV. INTERFERENCE IN FOCK STATES 

It is not evident that experimentalists can prepare a phase state as described above. It 
would seem more likely that they are working with Fock states, that is, states in which the 
particles numbers, A''^ and A^^, in the two condensates are known rather well. It seems likely, 
in any case, that one has more chance of initially preparing such a state. However, as several 
workers^"— have realized in recent years, and as we will show, that an interference pattern 
with some phase can still arise in a Fock state. Denoting the state sharp in particle number 
as lA*";!, Nh) and using the Bose annihilation relation 

Cr \Na,m) = sj^ (v^lAT, - l,Nb) + e^'^-'-v^lAT,, AT, - 1)) (24) 

the one-body density in a Fock state is 

D, = {Na,N,\clcr\Na,N,)=n (25) 

and of course there is no interference. Phase and particle number are conjugates, and the 
particle number is initially known. 

However, if we consider measuring the position of two particles simultaneously, some 
correlation should arise. We have 

Cr,Cr,\Na,N,) = ^ [^/ Na{Na - 1) \Na - 2, A^,) + e'''<'^+''^ y/ N,{N, - 1) {N,, m - 2) 

+ v'AyVb(e*('^-^i) + e*(''-'-^)) |(Ar, - l)(Ar, - 1))] (26) 

so that the particle number is now slightly less certain. The two-body Fock correlation 
function is 



D2 = {Na,Nb\cl/,^Cr,Cr,\Na,Nb) 

'KiK - 1) + N,iN, - 1) + KN, le^C^-'-i) + e*(^-^^ 



= — \nJN^ - n -^ NJN. - n -^ N^.N. |p*(k-ri) _^^i(k-r2)|2 

= n^[l-ha;cosk-(ri -rs)] (27) 

For two particles there is indeed a position correlation. As we make more and more measure- 
ments the state gets more and more mixed among states with various numbers of particles. 
We can rewrite Eq. ()27p in a somewhat different and useful way. A simple integration 
shows that 

D, = n' r^f[[l + xcos{k-r, + (j))]. (28) 



Compare this with Eq. ()22|) . We have a similar looking result except that we now integrate 
over all relative phases. Remarkably this result can be extended to higher order correlation 
functions. Indeed Ref. |lj shows that 

Dm = n"' ^ rr [1 + a; cos(k ■ r, + 0)] (29) 

Jo 2^t\ 

where it is assumed that m « N . 

In order to derive this result, we invert Eq. (jllj) . Multiply both sides by e"**^*^^"") = e~^'^^'' 

and integrate over to give 



Thus, by Eq. (fT^ . we find 

cr^ ■ ■ ■ cr, |iV„ N,) ~ / ^e-'^t>N, "Q ^(r„ 0) 10, N-m) (31) 

We will return to analyze this interesting state below. First consider the m-body Fock 
correlation function: 

D^ = {Na, Nt\ 4^ ■ ■ ■ C^Cr^ ■ ■ ■ Cr, |iV„ N^) (32) 

~ r ^ r ?e-(^-<^')^^ n A*(r„ 0')A(r„ 0) (0', iV - m| |0, iV - m) (33) 

Phase states are not actually orthogonal, but for large N they are essentially so as we show 
in Appendix A, so, if ?Ti -C A^, we can write (0', A^ — ?7i| |0, A^ — m) ~ 5(0 — 0') and 

f-2n 11 ™- 

D^- fX{\A{v,A)? (34) 

just as we claimed in Eq. (j29|) . We will show how to get this result by direct calculation, as 



done in Ref. |lj, in Appendix B. 

This result Eq. fj29|) has the form of Eq. ()23|) but with the unknown phase integrated 
over. That makes sense in that our initial Fock state did not have a phase defined and can, 
indeed, itself be expressed as a sum over phase states as in Eq. (jHUj) . 

Suppose we have started with a Fock state and have made m — 1 particle measurements 
and have found particles at positions Ri, ■ ■ ■ ,Rm-i and are about to measure the mth. 
Then the probability of finding the mth particle at position r^ is 

Pm = rr I ^9mi<f>) [1 + a; cos(k ■ r„ + 0)] (35) 
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where 

m—l 



9m{<p) = n [1 + ^ ^^^('^ ■ ^ + <^)] (^^) 



i=l 



As we show by direct simulation, g{(f)) develops a sharp peak at some a priori unpredictable 
phase. If one makes measurements from the first to the mth particle by this prescription, 
the peak becomes narrower as one proceeds. Of course, as more measurements are made, 
the particle number in each condensate becomes less certain (as, for example, in Eq. (j26j) ). 
so phase can be more sharply defined. 

Now look back at Eq. ()31|). After a fair number m of measurements, the real part of 
YYiLi ^(ti, 0) peaks up sharply at some (p value, call in 0q. That means that the measure- 
ments have converted the Fock wave function into a narrow sum of phase states around 
\(j)(^,N — m). The more measurements that are made, the better is the definition of the 
phase state. Measurements in a Fock state provide a way to prepare a phase state. One can 
understand the MIT experiments^ in this way. The starting state was prepared as two sepa- 
rate condensates, whose particle numbers could have been known; many subsequent particle 
measurements sharpened the phase to some random value and the final overall observation 
showed that phase. 

V. NUMERICAL SIMULATION 

We use Eqs. ()35|) and ()36p. One choses the initial position ri randomly, and then the 

next particle is chosen from the probability distribution P2 and so on. We will find that, if 

m is large enough, g^ peaks up at some a priori unpredictable phase angle, 0o which may 

fluctuate somewhat as m changes but gradually settles down. Starting a new experiment 

from the Fock state will lead to a randomly different phase angle. We will consider only 

the case in which the initial Fock state has Na = Nb = N/2, i.e., x = 1. It is convenient to 

Fourier expand Qm- That is, we write 

00 
5-^(0) = ao + ^ [ttg cos g0 + bg sin g0] (37) 

9=1 

In the integration of Eq. (jH3|) only oq, Oi^ and 61 will contribute and doing the integrals gives 

-Pm(rm) ~ 1 + -— cos(k ■ r„) - -— sin(k ■ r„,) (38) 

11 



If we define cos(0^) = ai/ ^Ja\ + h\, sin(0„) = hi/ ^a\ + h\, and A^ = \/ (A + &i/2ao, 
we can write 

Pm = -^'[1 + v4m(cosk- r™ cos (/)„,- sin k ■ r^sin*/)^)] = i^(l + y4„cos(k- r^ + 0^)) . (39) 

where i^ is a normalization factor, and 

tan(0^) = ^. (40) 

gives the the value of an angle in the mth experiment. Since P^ is a probability, we must 
have A^ < 1 so that P^ is always positive. Since A^ has that property we can write 
A^ = sin(am) where < a^ < t^ behaves like a polar angle. Then the emerging phase 
actually has a space angle designation (a^, 0„). We will find numerically that Am ~^ 1 
rapidly as we make measurements. In that case the probability of Eq. ()39|1 looks just like 
the density prediction of Eq. (jl} and moreover since gm{4') is a narrow function peaking at 
00, the phase defined by the Fourier coefficients is the same as that defined by the peak of 
gm, as seen using Eq. flH^ . 

We work in just one dimension. To choose from the probability of Eq. (J38|l . we need to find 
the cumulative probability Cm{x) = f^ dx'Pm{x') and then solve the equation R = Cm{x) 
for X, where P is a random number uniformly distributed on to 1. We take the box size 
L = 1 and choose a k value such that kL is an integer times 2tt to provide periodic boundary 
conditions. The normalization of the probability in Eq. ()39|) is just the factor K = 1/L = 1. 

Figs. 1 and 2 are plots of 0^ and A^ versus iteration number in a particular run of 200 
interations; each of these is found from the amplitudes oq, ai, bi at each iteration. There is no 
reason why Am should be unity from the outset. However, Am does always proceed to unity 
after a small number of iterations. The result is that 0^ approaches some sharply-defined 
random phase angle as predicted. Of course, for small m the fiuctuations are relatively 
large and settle down only after many measurements. This corresponds to an initially wide 
distribution, gm{<P), which, however, progressively narrows as more information is gathered. 
Fig. 3 plots the final angular distribution g2oo{4>)i it is indeed sharply peaked at the same 
value found from the iteration limit. Fig. 4 plots the final probability distribution, Eq. ()39j) . 
versus postion x and also shows a histogram of the positions found in the 200 iterations on 
the same scale. One sees that these positions really do fall in the given distribution with 
the expected oscillations and with the same phase as gotten in the two other ways. 
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VI. DISCUSSION. 

We have shown how one might rigorously treat the interference of two Bose condensates. 
The usual assumption of two condensates with individually known phases gets involved with 
the thorny questions about whether one can usefully define the phase of a single condensate. 
However, using such states does lead directly to the usual fomulas for interference patterns 
based on very simple assumptions. This procedure remains unsatisfying since it may not be 
very obvious to the experimenter how to prepare such phase states before looking at the in- 
terference. Experimentally it seems to make no difference, since without special preparation 
experiments, even with Fock states, we have seen how an interference pattern arises using 
Bose condensates. The discussion of Sec. IV shows explicitly why such preparation was not 
necessary. Even if one starts with a state where the particle number in each cloud is precisely 
known and many particles are involved in the measurement, one finds a perfect interference 
pattern emerging, with a well-defined relative phase. Starting from a state having a definite 
number of particles, the experiment will end up with a state having a quite definite value of 
the relative phase. Thus this procedure actually provides a method of preparing the phase 
state discussed in Sec. III. Starting from a Fock state, make, say, two hundred position 
measurements to get to a narrow g2oo{4')] the phase of the wave function of the remaining 
state of A^ — 100 total particles will be well-defined. The final result is likely a phase state 
with a known total number of particles such as that discussed at Eq. (Illj). but an unknown 
number in each condensate. 

In theoretical superfluid calculations it is simpler to treat the problem with definite phases 
than to use a Fock state. However, the actual existence of such a broken symmetry state 
is subject to question.—"— In a ferromagnet one actually has broken the symmetry of the 
various directions of magnetization from the presence of some small external field. The 
existence of a field that would make (ip) non-zero is not so clear, since such states do 
not conserve particle number. The existence of a well-defined relative phase established by 
measuring particle positions can be established without worrying about broken symmetry. 

If a person feels uncomfortable with the idea of a phase emerging from the series of 
measurements done on particle position, then he or she might assume, with no change in 
theoretical prediction, that that phase pre-existed within the clouds of particles. That is, 
the individual condensates had some relative phase, before they met, as a so-called "hidden 
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variable" and the experiments simply bring out this previously hidden phase. Of course, 
in the next realization of the experiment, starting again from a Fock state, the phase will 
surely emerge with a randomly different value, in accordance with conventional quantum 
mechanics, which expresses the Fock state as a sum over all phase states as given in Eq. 

dnoi). 

APPENDIX A. COHERENT STATES 



Consider the following wave function to describe a single momentum state k, which is 
made up of a mixture of states of known particle number N^. 



I«fc) 



e 2 



hi 






a 



Nk 



fNl 



\Nu). 



(41) 



This function is properly normalized. The parameter a^ is complex and we separate it into 
magnitude 7^, and phase 0^ according to the notation 



«fc = 7fce 



(42) 



We can easily compute the average number of particles A''^ in this state. Let at be the 
destruction operator for particles in state \N].) . Then 



2JVfc 



A^, = {a,\ aW W) = e-^l ^ ^AT, = 7^ 



(43) 



Nk 



and 7fc = \ak\ = V^ • 

The state \ak) has the nice property that it is an eigenstate of the lowering operator a^ : 



Ofc \ak) 



< 



Nt: 



J2^ak\Nk) = ake--lYl 



a 



Nk-l 



Nk 



m. 



Nt 



ViNk - 1) 



\Nk -l) = ak \ak) . (44) 



Thus Ofc has a non-zero expectation value in this state: 

{ak\ak\ak) = ak = ^/Wke"^'^ (45) 

Clearly the states \ak) provide a definite phase. Of course they are then not eigenstates of 
the number operator. 

Next construct a multi-level many-body state with many k values possible. This takes 
the form la^^, ^^2 a^a " " " ) = l«fci) \0ik2) bfcg) ■ ■ ■ 



\{ak}) = ittfci, ttfca «fc3 ■ ■ • ) = x^ n 

{Nk} k 



e 2 



H^^ 



VW. 



\Nk,,Nk,,Nk,---) 



(46) 
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where {Nk} means sum over all possible numbers of particles A^^. in all the /c-states. 

With such a state, we can consider the expectation value of the full field operator -^(r) 
Expand the field operator in plane wave states as 



^w = \/^E^^''^«^ (47) 



where V is the volume of the system. We have 

({a.}| ^(r) \{a,}) = J^ ^J^e^'^e^^-. (48) 

k 

If one of the fc-states is macroscopically occupied, say, the momentum state k = ka, then we 
can write 

'-^{r) \ = y^e^^^-^'e''^^ + ^ (49) 

where Ha = Na/V and ^ is the total contribution of the non-condensed states. The leading 
term '^/'^(r) = ^/n^e^^°■'^e^^'' represents a condensate wave function having a definite phase 
(/)„ but non-sharp number of particles. 

Consider the case of the interference of a double condensate in momentum states ka and 
kh. With coherent states having only these two momentum states occupied we can write 

Ka,) = e-K^^+^i) 5^ ^1^ \Na, N,) (50) 

where the averages Na = 7^ and Nf, = 7^ are both macroscopic quantities. We manipulate 
the sums slightly in terms of particle creation operators a^ and b^ for the two states. If 
N = Na + Nfy then we have 

K«.) = e-K^^+^?) Y: /";""""" ^\Na,N-Na) 
^^^NJiN-NaV. 

N,Na "■ ^ ^^' 



= e-K^^+7n^_l(«,«t + «^5t)^|0) 

TV 

The second last form of Eq. (0) is the "phase state" used in Sec. IV. 

\aaab; N) ~ {aaa^ + abb'^)^ |0) . (52) 
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We see that it is a substate of the more general coherent state. 



APPENDIX B. NEAR-ORTHOGONALITY OF PHASE STATES 



We calculate the inner product of two phase states to show that they are nearly orthogonal 
for large particle number. From Eq. (jll|) we find 



N 



giM z_^ n\[N — ny. 
This is a very sharply peaked function of (0 



^2gi(<^-</,') 



1 N-n 



1 



1 + 7 e' 



2i(4>-<t>') 



iN 



(53) 



b') . To see this, Taylor expand the logarithm 



of this in powers of 
. The result is 



b') and then exponentiate the result keeping only terms to 



b',N\(j),N) =exp 



A^ 



exp 



N 



[l + l' 



2(1 + ^2)2 V^ ^ 

In the limit of very large A^, this is proportional to a delta function of the 
assumed in the discussion of Sec. IV. 



(54) 

6') as we 



APPENDIX C. ALTERNATIVE DERIVATION OF THE D^ EQUATION 

We want to derive the general expression of Eq. (j29|) for the correlation function D^ 
Consider this quantity in its original form for a Fock state: 



Dr 



{Na,N,\cl^---cl 



■Cr,\Na,Nh) 



— {Na, Nt,\ (at + e-^'^-'-^fet) ...(«! + e-i^-^fot) 



X 



(a + e'^-^-^b) • ■ ■ (a + e'^'^'b) \Na, N^,) 



(55) 



Because this is a diagonal matrix element in Fock space, each time an a occurs, there must 
be a matching a^. Similarly for the b operators. We are assuming m <^ Na oi Nb so that 
we can always write a\Na — p,Nh — I) ~ y/N^ \Na — p — l,Nh — I) , etc. Thus each a'^a gives 
Na, and each b'^b an Nt. Consider a particular combination product: 



(a^ + e-'^''^b^) {a + e'^'^'b) = a^a + b^b + a^be^^''' + b^ae 



be ^''■^^ 



(56) 



with the restriction that every time an e*'''''-type term occurs there must be a cor- 
responding e~*'^'"^-type term somewhere in the overall product to give the proper bal- 
ance of creation and destruction operators. Thus one gets a series of terms of the form 
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Fq^{rm)Fq^_,{rm-i) " " " F,^ (rs)^^, (fi) where 

Fo(r,) = Na + Nt 



F±i{r,) = ^NJ^e^^'^-'^ (57) 

and the sum of all the qj vanishes. That is, we have 

where {g} means sum on all qi with the restriction that ^^ qi = 0. 
The restriction on the q values can be lifted if we insert the integral 

f'S^«'-^E,„o (59) 

which allows us to write 

= rT ^ n [1 + xcos (k ■ r, + 0)] (60) 

Jo SttA^A 

as we wished to prove. 
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FIG. 1: The phase angle (p^ as a function of interation step. 
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FIG. 2: The amplitude A = sin a as a function of interation step. This always proceeds to unity. 
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FIG. 3: The angular distribution g{(p) as a function of angle. This peaks at the same phase angle 
as given in Fig. 1. 
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